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1.  INTRODUCTION 


Reported  herein  is  a summary  of  research  on  elastic-plastic  analysis 
of  cracks,  during  the  period  February  1,  1974  to  February  1,  1978,  under 
the  support  of  AFOSR  Grant  74-2667. 

The  research  described  has  been  performed  at  the  School  of  Engineering 
Science  and  Mechanics,  Georgia  Institute  of  Technology,  under  the  direction 
of  Satya  N.  Atluri.  Three  doctoral  theses  were  written,  by  Drs.  Nakagaki, 
Kathiresan,  and  Chen,  based  on  research  performed  under  this  grant;  these 
are  referenced  at  the  end  of  this  report.  Dr.  Nakagaki  later  participated 
in  this  research  as  a Post-Doctoral  Fellow  and  made  significant  contributions. 

Dr,  W.  J.  Walker  of  AFOSR  provided  invaluable  encouragement  and  has 
aided  immensely  in  defining  the  practical  objectives  of  the  basic  research. 

2.  AREAS  OF  RESEARCH 

The  studies  under  the  present  grant  included:  (a)  development  of  two- 
dimensional  hybrid  finite  element  procedures  to  calculate  stress  intensity 
factors  corresponding  to  general  r stress  singularities  in  isotropic 
as  well  as  anistropic  materials;  (b)  study  of  the  convergence  of  the 
assumed  displacement  hybrid  finite  element  procedure  in  fracture  mechanics 
problems;  (c)  development  of  a three-dimensional  hybrid  finite  element 
procedure  to  calculate  the  elastic  combined  mode  stress-intensity  factors 
Kf,  ^III  along  an  arbitrarily  curved  three-dimensional 

crack  front;  (d)  hybrid  finite  element  solutions  of  fundamental  three- 
dimensional  crack  problems;  (e)  development  of  a two-dimensional  finite 
element  procedure  for  analyzing  plane  problems  of  fracture  involving  large- 
scale  yielding  conditions  and  under  cyclic  loading;  (f)  study  of  the  J- 
integral  as  a ductile  fracture  initiation  condition;  (g)  studies  of  finite 
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deformation  effects  near  the  crack-tip;  (h)  analysis  of  stable  crack  growth 
under  rising  load  and  study  of  criteria  for  loss  of  stability  of  growth  in 
ductile  materials;  and  (i)  elastic-plastic  analysis  of  effects  of  crack 
closure  on  fatigue  crack  growth  rates. 

A brief  report  on  the  salient  results  obtained  in  each  of  the  above 
studies;  a listing  of  publications  that  have  resulted  from  this  research; 
papers  presented  at  professional  conferences;  and  seminars  given  at  other 
organizations,  are  given  in  the  following. 

3.  SUMMARY  OF  RESEARCH  ACCOMPLISHMENTS 
u Two-Dimensional  Linear  Elastic  Fracture  Analysis. 

A basic  procedure,  based  on  an  assumed  displacement  hybrid  finite 
element  method,  to  calculate  the  stress  intensity  factors  for  mixed  mode 
behaviour  of  arbitrary  shaped  curved  cracks  in  plane  stress  and  plane 
strain  problems  involving  anisotropic,  nonhomogeneous , but  linearly  elastic 
materials  has  been  developed  [Refs.  1,  4].  The  attendant  finite  element 
program  has  been  used  to  study:  (i)  unsymmetric  boundary  conditions  and 
configurations  such  as  oblique  edge  cracks  in  tension  plates;  (ii)  general 
plane  fracture  problems  with  stress,  displacement,  or  mixed  conditions; 
(iii)  curved  cracks  through  the  thickness  of  tension  specimens;  (iv) 
problems  involving  material  inhomogeneity  such  as  cracks  intersecting  bi- 
material interfaces  in  tension  plates,  wherein  stress  singularities  of 
- a 

the  general  r type  may  arise;  (v)  cracks  meeting  a free  surface  such  as 
single-edge-notched  tension  plates  with  cracks  of  varying  aspect  ratio  to 
study  the  effects  of  free  surfaces  on  crack-tip  stress  intensity  factors; 
(vi)  cracks  in  orthotropic  sheets,  with  the  crack  line  located  either 
symmetrically  or  unsymmetrically  with  respect  to  directions  of  principal 
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£ixes  of  orthotropy,  and  for  various  values  of  Ex/Ey  to  obtain  stress 
intensity  factors  for  laboratory  test  specimens.  These  results  are  document- 
ed in  Refs.  5-8. 


B.  Study  of  Convergence  of  the  Hybrid  Element  Procedure  for  Crack  Problems. 

A conyergence  study  was  conducted  for  the  assumed  displacement  hybrid 
finite  element  method  for  solving  stress  intensity  factors  in  linear  fracture 
mechanics.  Also  methods  to  improve  the  satisfaction  of  stress-free 
conditions  on  a boundary  segment  of  an  element  adjoining  the  crack  surface 
were  developed.  It  was  found  that  a proper  way  to  obtain  convergent 
solutions  for  K-factors  is  to  keep  the  size  of  the  hybrid  crack  elements 
near  the  crack-tip  at  constant,  pre-determined  optimum  size,  while  the  size 
of  the  surrounding  regular  elements  is  progressively  decreased.  These 
results  are  documented  in  Refs.  1,  5,  6 and  7. 

C.  Three-Dimensional  Hybrid  Crack  Element  for  Linear  Elastic  Fracture 

Analysis.  , 

1 

A finite  element  procedure  was  developed  for  the  calculation  of  • 

modes  1,  II  and  III  stress  intensity  factors,  which  vary  along  an  arbitrarily 
curved  three-dimensional  crack  front  in  a structural  component.  The  present 
finite  element  model  is  based  on  a modified  variational  principle  of  potential 
energy  with  relaxed  continuity  requirements  for  displacements  at  the  inter-  ^ 

element  boundary.  The  variational  principle  is  a three-field  principle, 
with  the  arbitrary  interior  displacements  for  the  elements,  interelement 
boundary  displacements,  and  element  boundary  tractions  as  variables.  The 
unknowns  in  the  final  algebraic  system  of  equations,  in  the  present  model, 
are  the  nodal  displacements  and  the  three  elastic  stress  intensity  factors. 

Special  elements,  which  contain  proper  square  root  and  inverse  square  root 
variations  in  displacements  and  stresses,  respectively,  are  used  in  a fixed 
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region  near  the  crack  front.  Interelement  displacement  compatibility  is 
satisfied,  by  assuming  an  independent  interelement  boundary  displacement 
field,  and  using  a Lagrangian  multiplier  technique  to  enforce  such  inter- 
element compatibility.  These  Lagrangian  multipliers,  which  are  physically 

the  boundary  tractions,  are  assiuned  from  an  equilibrated  stress  field  de-  j 

rived  from  three-dimensional  Beltrami  (or  Maxwell-Morera)  stress  functions 

that  are  complete.  However,  considerable  care  had  to  be  exercised  in  the 

use  of  these  stress  functions  such  that  the  stresses  produced  by  any  of 

these  function  components  are  not  linearly  dependent. 

Since  the  method  is  based  on  a rigorous  variational  principle,  which 
enforces  at  least  on  an  average  the  conditions  of  interelement  displacement 
continuity  when  v'^  type  displacements  are  included  in  the  near-tip  region, 
the  convergence  of  the  finite  element  solution  for  nodal  displacements  as 
well  as  the  stress  intensity  factors  is  established  mathematically. 

The  geometry  of  the  "basic  element"  used  presently  is  a 20-node 
"isoparametric"  brick  element,  with  60  degrees  of  freedom  per  clement.  i 

Two  options  were  developed:  (a)  a singular  element  for  use  near  the  crack 
front,  wherein  the  stress  Intensity  factors  K^,  and  are  constant, 

and  (b)  a singular  element  where  in  the  intensity  factors  Kj,  and 

themselves  vary  quadr iat ical ly.  It  is  believed  that  this  option  will 
reduce  the  number  of  elements  to  be  used  near  the  crack  front,  in  order 
to  obtain  an  accurate  variation  of  the  intensity  factors  along  the  crack 
front . 

The  relevant  matrices  were  evaluated  numerically,  using  non-product  | 

type  quadrature  formulae,  with  proper  mathematical  transformations  being  i 

I 

used  when  singular  type  functions  arc  encountered  in  stresses  and  strains  j 

in  the  near-tip  region. 
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The  mathemat ical  development  of  the  3-D  hybrid  crack  element  and  related 
computational  details  are  documented  in  Refs.  2,  10  and  13. 

D.  Hybrid  Finite  Element  Solutions  of  Fundamental  Three-Dimensional  Crack 
Problems . 

The  basic  development  for  3-D  hybrid  crack  element  discussed  in  (C)  above 
has  been  exhaustively  tested  in  analyzing  various  basic  problems  of  3-dimen- 
sional  fracture  that  are  of  current  Interest.  These  include:  (i)  semi- 
elliptical  surface  flaws  of  various  aspect  ratios  and  various  crack-depth 
ratios  in  plates  under  tension  as  well  as  bending;  (ii)  semi  elliptical 
surface  flaws  of  various  aspect  ratios  and  various  depth  ratios  on  the  inner 
and  outer  surface  of  pressurized  thick  shells;  (iii)  semi-elliptical  surface 
and  quarter-elliptical  corner  cracks  near  fastener  holes  in  plates,  typical 
of  aircraft  structural  configurations.  Moreover,  the  method  was  completely 
checked  as  to  its  accuracy  and  convergence  in  several  other  problems  such 
as  Sneddon's  embedded  penny  shaped  crack,  buried  elliptical  cracks,  and 
several  other  problems,  for  which  analytical  solutions  exist.  Also,  the 
question  of  the  optimum  size  of  a singular  element  near  the  crack  front  has 
been  studied  thoroughly.  All  these  results,  along  with  the  detailed  mathe- 
matical derivations,  have  been  presented  in  a recent  Ph.  D.  thesis  by 
Kathiresan  [2],  and  are  submitted  to  AFOSR  for  approval  for  publication  as 
an  AFOSR  Technical  Report  [21].  The  results  for  various  problems  have  also 
been  published  in  open  literature  [13,  17,  19,  24,  26,  28,  31,  32]. 

More  recently,  the  3-D  hybrid  crack  element  procedure  has  been  applied 
to  study  the  crack  growth  behaviour  in  adhesively  bonded  structure,  wherein, 
each  layer  was  treated  in  a three-dimensional  fashion.  Specifically  results 
were  obtained  for  a crack  through  the  thickness  of  only  the  outer  layer,  near 


a hole  through  the  thickness  of  the  entire  stack  of  layers.  These  are  believed 
to  be  the  only  fully  three-dimensional  analysis  results  for  the  problem  under 
consideration;  and  are  presented  in  Ref.  [32]  along  with  a discussion  as  to 
the  validity  of  the  obtained  results. 

K . Two-Dimensional  Hybrid  Element  Procedure  for  E las t ic -Plas t ic  Analysis  of 
Large  Scale  Yielding  Fracture: 

To  study  the  problem  of  ductile  fracture  under  large-scale  yielding 
conditions,  research  was  conducted  in  formulating  and  developing  a second 
generation  finite  element  procedure  with  the  objectives;  (i)  developing 
circular-sector  shaped  embedded -s ingu lar ity  rinite  elements  near  the  crack- 
tip.  The  correct  r dependence  of  the  diminant  singular  solution,  correspond- 
ing to  the  nonlinear  material  model  (Ramberg-Osgood  law),  was  embedded  in 
these  near  tip  elements,  where  as  the  ©-dependence  was  approximated  in  each 
sector  element  and  solved  for  in  the  sense  of  the  finite  element  method; 

(li)  maintaining  continuity  of  displacements  and  tractions,  between  near-tip 
elements  with  singular  stress/strain  assumptions  and  the  far  field  elements 
with  regular  stress/strain  assumptions,  through  a hybrid  displacement  finite 
element  model,  (iii)  using  a J2-flow  theory  of  plasticity  and  arbitrary 
kinematic  hardening  wlilch  will  accurately  model  the  Bauschinger  effect  under 
fully  reversed  and  cyclic  loading,  (iv)  using  an  incremental  finite  element 
solution  procedure  that  will  be  suitable  in  the  limiting  case  of  elastic- 
perfect-plastic  materials,  as  well  as  cyclic  loading  situations.  The  "initial- 
stress"  iteration  approach  was  used  for  this  purpose,  (v)  developing  a more 
accurate  finite  element  method  for  incremental  analysis  of  elastic-plastic 
problems.  The  more  common  approach  in  leterature  is  to  use  "constant ■ stress" 
elements,  and  based  on  stress  level  in  each  element,  t!ie  wliole  element  either 
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yields  or  stays  elastic.  Thus  in  problems  such  as  the  present,  where  the 
yield  zones  near  the  crack  tip  play  a dominant  role  in  the  analysis  and 
its  interpretation,  in  order  to  obtain  a reasonably  accurate  description  of 
the  yield  zone,  a very  fine  finite  element  mesh  is  needed.  However,  if 
higher  order  elements  are  used,  and  if  "plasticity-correction"  iterations 
are  performed  at  several  points  within  the  element,  it  then  becomes  possible 
to  give  a smoother  definition  of  the  yield  zone.  Thus,  a portion  of  the 
element,  in  the  present  formulation,  can  yield  while  the  rest  of  the  element 
can  remain  elastic. 

To  start  with,  crack-tip  blunting  and  other  finite  geometry  changes 
near  the  crack-tip  were  ignored,  and  thus  the  analysis  was  restricted  to  a 
small-deformation  assumption.  This  procedure  was  used  to  solve  a hypothetical 
3-point  bend  test  specimen  which  was  suggested,  as  a standard  test  for  analyti- 
cal solutions,  by  the  ASTM  Committee  E-24  Round  Robin  Group  on  Elasto-Plastic 
Fracture  Criteria.  A detailed  analysis  of  these  results,  with  various  fracture 
criteria  in  mind,  was  presented  in  (Atluri  and  Nakagaki  [ 12 j , and  Atluri  and 
Nakagaki  [14j.  However,  no  other  independent  solution,  from  members  of  the 
Round  Robin,  with  which  to  compare  our  results  is  available  at  the  time  of 
this  writing. 


F.  Study  of  the  J-integral  as  a Fracture  Initiation  Criterion  Under  Large- 
Scale  Yielding. 

The  above  small-deformation  analysis  procedure  was  used  to  analyze  a 
compact  tension  cracked  specimen,  of  A533B  steel,  for  which  experimental 
fracture  data  was  reported  by  the  Westinghouse  Research  Laboratories.  The 
thickness  of  the  experimental  specimen  was  roughly  1/2-1/3  of  the  characteristic 
inplane  dimensions.  Thus  it  was  not  a priori  certain  whether  the  conditions 
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near  the  crack-lip  can  be  char. -c t er i zed  mathematically  as  those  of  plane  stress 
or  of  plane  strain.  'I'luis  the  prolilem  was  analyzed  using  two  types  of  two- 
dimensional  approximations.  Also,  tlie  experimental  uniaxial  s t ress /s l ra in 
curve  supplied  to  us  by  the  Wes t Ingliouse  Laboratories  was  seen  to  possess  a 
yield  point  inst  .-ibi  1 ity , whlcli  the  numerical  procedure  could  not  account  for. 
riuis , two  different,  smootli,  K.aniberg  Osgood  type  curves  were  used  to  fit  I l>e 
actual  experimental  stress-strain  curve.  l•'rom  tliese  analyses,  it  was  found 
th.at:  ( i ')  the  mathematical  characterization  of  the  ma  t er  ia  I -proper  t y data 

had  no  significant  effect  on  the  results,  (il^  the  empirical  formula  of 
Rice  et  al.  for  estimating  ,t  was  in  significant  error  in  the  case  of  llte 
compact  tension  specimen,  and  (ill)  the  computed  nl.ane-stress  values  for  .1 
were  muclt  closer  to  the  experimental  data,  as  compared  to  the  plane -strain 
values  wliicli  differed  by  a factor  of  about  3.  Thus  it  was  concluded  that 
the  plane  stress  conditions  more  accurately  cliaracter  i ze  llie  plastic  flow 
near  tlie  crack  tip  in  "small"  specimens.  Tliese  results  are  documented  in 
detail  in  (Atluri  and  Nakagaki  [ ^ • 

However,  even  the  computed  plane  stress  values  for  .1,  based  on  the 
present  analysis  with  a "non-growing  stationary  crack"  assumption,  differed 
by  about  13-20'/.  from  the  cited  experimental  results  for  .1  at  fracture,  from 
subsequent  private  communications  with  Regley  and  l.andes,  who  conducted  the 
experiments,  it  was  learned  that,  in  fact,  there  may  have  been  stable  crack 
growth,  which  was  not  monitored,  in  the  experiments. 

Taking  into  account  the  effect  of  this  stable  crack  growth,  and  the 
results  of  a refined  and  consistent  finite  deformation  analysis  of  the  problem, 
as  discussed  in  ((',)  below,  it  was  concluded  that  ,1  is  a computationally 
advantageous  criterion  for  the  onset  of  crack  growth  in  ductile  materials. 
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G.  Studies  of  Finite  Deformation  Effects  near  the  Crack-tip. 


To  delineate  the  effects  of  crack-tip  blunting,  and  the  attendant  finite 
geometry  changes  near  the  crack-tip,  under  large  scale  plastic  yielding 
conditions,  a finite  deformation,  embedded  singularity,  elastic  plastic 
Incremental  finite  element  procedure  was  developed  [16j.  The  test  problem 
was  that  of  a 3-point  bend  fracture  test  specimen,  the  experimental  data 
for  which  was  reported  and  made  available  to  us  by  Westingliouse , A detailed 
comparison  of  the  results,  for  stresses,  strains,  displacements,  and  J,  from 
tills  finite  deformation  analysis,  and  those  obtained  using  the  above  small- 
deformation  analysis  was  made.  Significant  conclusions  of  the  above  compari- 
son were:  (i)  not  with-standing  the  crack-tip  blunting,  J still  remains 
as  a valid  parameter  to  characterize  the  severity  of  the  near-tip  condition, 

(ii)  since  the  well  known  Rice's  definition  of  J is  invalid  for  finite- 
deformation,  a modified  definition  must  be  derived.  This  new  definition 
involves  the  strain  energy  density  function  that  is  dependent  on  Piola- 
Lagrange  (unsymmetric)  stress  tensor,  and  the  traction  term  appearing  in  the 
.1- Integral  should  be  interpreted  as  Piola  Lagrange  tractions  in  the  deformed 
geometry  as  referred  to  the  undeformed  configuration.  Wlien  this  modified 
definition  for  J is  employed,  a more  accurate  path-dependence  (within  + 1.5%') 
was  noted  for  J,  as  compared  to  the  smal 1 -deformat  ion  case  wherein  the  J 
value  on  paths  closest  to  the  crack-tip  was  about  15%  lower  than  that  on 
paths  in  the  far  field.  A detailed  discussion  of  these  results  was  presented 
in  (Atluri,  Nakagaki,  and  Chen  [16]. 

One  interesting  observation  was  that  at  all  load-point  displacement  levels, 
including  that  at  fracture,  there  was  an  excellent  correlation  (within  + 47.) 
between  the  presently  computed  results  for  J and  those  in  the  Westlnghouse 
experiments.  This  suggested  that  no  appreciable  stable  crack  growtlt  may 
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have  been  present  in  the  experiment.  Subsequent  private  communications  with 
Dr.  Landes  at  Westinghouse  appeared  to  confirm  this.  Thus  because  of  tlie 
particular  specimen  geometry  and  loading  conditions,  3-point  bend  specimen 
may  offer  the  advantage  of  precisely  measuring  for  crack  growth  initiation 

for  a given  ductile  material. 

Further  refinements,  such  as  the  incorporation  of  a "knee-correction", 
to  achieve  better  convergence  of  plastically  adjusted  stresses,  and  a more 
accurate  incremental  elastic-plastic  law,  were  made  in  the  above  finite 
deformation  analysis.  The  details  of  the  nature  and  form  of  singularity  of 
actually  computed  strains  and  stresses,  and  their  variation  in  the  angular 
coordinate  near  the  crack  tip  were  critically  reviewed  and  compared  with 
the  small-deformation  analysis,  using  deformation  plasticity  theory,  of 
Hutchinson,  Rice  and  Rosengren.  These  results  are  discussed  in  detail  in 
(Atluri  and  Nakagaki  [20]).  Because  of  the  finite  root  radius  of  the 
blunted  notch,  there  can  only  be  a strain  concentration,  however  large,  at 
the  tip  of  the  notch.  The  effect  of  this  strain  concentration,  as  opposed 
to  the  strain -singularity  at  the  blunted  crack  tip,  was  properly  accounted 
for  in  the  analysis  of  [20]. 

The  refined  finite  deformation  analysis  procedure  was  also  ued  to 
analyze  three  fracture  test  specimens,  of  3-point  bend,  center -cracked , and 
compact  tension  types,  of  configurations  identical  to  those  in  the  experiments 
of  Westinghouse  Laboratories,  The  obtained  numerical  data  was  analyzed  with 
both  the  J and  COD  criteria  in  mind,  and  the  following  conclusions  were 
drawn:  (i)  the  computed  J correlated  well  with  the  cited  experimental 

results  for  all  the  specimens;  (ii)  J was  found  to  be  directly  correlated 
with  COD  for  these  test  specimens.  The  above  data  indicated  that,  for 
these  specimens,  J = 1.44o'  (COD)  for  moderately  hardening  materials,  where 
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is  the  yield  stress.  Thus  it  can  be  seen  that  (X)D  may  also  be  used  as  a 
fracture  criterion;  (iii)  liowever,  it  was  found  that  the  computation  of  COD, 
as  essentially  a near-tip  geometric  quantity,  for  arbitrary  plane  problems 
of  ductile  fracture,  is  not  well  defined.  \ttempts  to  correlate  COD  with 
sucli  near  tip  geometrical  quantities  such  as  the  crack  opening  at  points 
where  the  elastic-plastic  boundary  intersects  the  crack  profile  were  found 
to  be  discouraging;  (iv)  the  empirical  formula  of  Rice  et  al.  for  J was  once 
again  found  to  be  inaccurate  for  the  compact  tension  specimens.  These 
results  are  elaborated  upon  in  (Atluri,  Nakagaki,  and  Clien  [18]). 

Once  again,  the  computed  J for  the  compact  tension  specimen. based  on  the 
finite  element  modeling  of  a non-growing  crack,  were  found  to  be  about  10-157„ 
less  than  in  the  experiments.  Again  this  suggested  that  there  may,  in  fact, 
have  been  a subcritical  crack  growth  prior  to  fracturing  in  the  experiments, 
thus  providing  motivation  for  our  study  (H)  on  stable  crack  growth,  described 
below. 

It . Analysis  of  Stable  Crack-Growth  and  Stiidv  of  Criteria  Governinit  Loss  of 
Stability. 

A finite  element  methodology  was  developed  to  study  tlie  phenomena  of 
stable  crack  growth  in  two-dimensional  problems  involving  ductile  materials. 

Crack  growtli  is  simulated  by  (it  the  translation  in  steps,  of  a core  of  sector 
elements,  with  embedded  singularities  of  Hutchinson-Rice -Rosengran  type  by  an 
arbitrary  amount,  Aa  in  each  step,  in  the  desired  direction;  (ii)  reinter- 
polation of  the  requisite  data  in  the  new  finite  element  mesh;  and  (iii) 
incremental  relaxation  of  tractions  in  order  to  create  a new  crack  face  of 
lengtli  Aa.  Steps  (i)  and  (ii)  were  followed  by  corrective  equilibrium- 
check  iterations.  A finite  deformation  analysis,  based  on  the  incremental 
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updated  Lagrangean  formulation  of  the  hybrid-displacement  finite  element 
method,  is  used. 

The  present  procedure  is  used  to  simulate  available  experimental  data  on 

stable  crack  growth,  and  thus  study  the  variation,  during  crack  growtl\,  of 

certain  physical  parameters  that  may  govern  the  stability  of  such  growth 

and  the  subsequent  onset  of  rapid  fracture.  Attention  is  focussed  in  this 

/Va 

study  on  the  parameters;  G the  energy  release  to  the  crack-tip  per  unit 
crack  growth,  for  growth  in  finite  steps  Aa.  calculated  from  global  energy 
balance  considerations;  the  energy  release  to  a finite  "process  zone"  near 

the  crack-tip  per  unit  crack  growth,  for  growth  in  finite  steps  Aa,  calculated 
again  from  global  energy  balance  considerations;  and  the  crack  opening  angles. 

However,  only  the  first  phase  of  our  study,  viz.,  the  simulation  of 
available  experimental  data,  has  been  completed. 

It  is  recognized  that  formulating  any  criterion  or  criteria,  governing 
the  loss  of  stability  of  crack-growth,  based  on  numerical  simulation  of  a 
few  experimental  data  is,  at  best,  a risky  proposition.  Thus,  we  defer  any 
conclusions  regarding  such  criteria  until  the  completion  of  second  phase  of 
our  researcl),  to  be  conducted.  However,  the  results  obtained  so  far.  can 
lead  to  the  following  conclusions  that  may  be  germane  to  the  problem  of  stable 
crack  growth  in  ductile  materials:  1)  A direct  numerical  proof  is  provided 

* A 

for  the  original  hypothesis  of  Rice  that  G -•  o as  Aa  — o for  those  materials 
for  which  the  flow  stress  saturates  at  a finite  value  of  large  strain.  Thus, 
for  any  meaningful  numerical  study  of  stable  crack  growth,  a finite  growth 
step  must  be  postulated.  Results  obtained  during  this  research  and  reported 
in  [27,30]  may  be  useful  in  providing  guidelines  for  choosing  Aa  such  that 

A 

the  numerically  computed  G is  not  sensitive  to  tlie  errors  inherent  in 
numerical  processes  such  as  in  the  finite  element  method.  2)  Since  tlie 
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magnitudes  of  G “ and  Gp  are  clearly  shown  [27,  30j  to  depend  on  the  post- 
ulated magnitudes  of  "finite"  growth  steps  Aa,  in  the  finite  element  model- 
ing, it  is  clear  that  any  criteria  governing  the  loss  of  stability  of  crack 
growth  cannot  be  based  on  the  absolute  magnitudes  of  these  quantities.  Any 
such  criteria  can  only  be  based  on  the  relative  qualitative  behaviour  of 

*A  -Va 

G and  Gp  for  the  postulated  growth  step  Aa.  3)  In  view  of  the  previous 
observation  and  the  fact  G and  Gp  vary  substantially  during  the  crack 
extension  process,  it  is  clear  that  the  generalizations  of  Griffith's  approach, 
in  tlte  sense  that  a ductile  material  has  some  characteristic  work  of  separa- 
tion, per  unit  new  crack  area,  and  that  this  is  to  be  equated  at  the  critical 
condition  to  the  rate  of  surplus  of  work  done  on  tlie  material,  cannot  be  made 
in  situations  of  stable  growth  under  large  scale  yield  conditions.  41  The 

only  discernible  trends  near  the  points  of  fracture  as  observed  in  experiment, 

*A 

for  both  the  cases  studied,  is  the  marked  change  in  the  behaviour  of  G , 

*A 

CTOA,  and  to  an  extent  in  Gp  , near  these  points,  when  tliese  quantities 
reverse  their  nomotonical ly  increasing  trend  during  the  prior  extension 
process.  Wliile  a theoretical  argument  explaining  tills  is  lacking,  it  yet 
remains  to  be  seen  whether  these  observations  can  be  used  to  numerically 
predict  loss  of  stability  of  growth  in  different  specimens  of  the  same 
material.  This  is  the  object  of  our  work  yet  to  be  conducted. 

A detailed  discussion  of  the  present  analysis  of  stable  crack  growth 
and  results  leading  to  the  above  conclusions,  have  been  presented  in  (Atluri 
and  Nakagaki  [27],  and  Nakagaki,  Oien,  and  Atluri  [30]). 

I . Analysis  of  Effects  of  Crack-Closure  on  Fatigue  Crack  Growth  Rates. 

The  crack  growth  simulation  procedure  described  in  (G)  above  was  also 
used  to  study  the  crack-closure  and  opening  stresses  under  Mode  1 type 
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constant  amplitude  fatigue  loading.  The  values  of  closure  and  opening  stresses 
for  a center  cracked  Aluminum  panel  subject  to  constant  amphitude  cyclic 
loading,  which  causes  a "small-scale"  yielding  near  the  crack-tip,  were  found 
to  correlate  well  with  existing  results  in  literature  for  Identical  problems. 
However,  the  number  of  finite  element  degrees  of  freedom  used  in  the  present 
study  was  found  to  be  an  order  of  magnitude  lower  than  in  current  literature 
for  comparable  accuracy.  These  results  are  yet  to  be  published  [30]. 

Considerations  of  general  spectrum  loading  to  understand  the  crack 
growth  acceleration  and  retardation  effects  still  remains  and  is  the  object 
of  our  research  yet  to  be  conducted. 
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Chen,  "On  Hybrid  Finite  Element  Models  for  Linear  and  Nonlinear  Fracture 
Analyses",  in  Numerical  Methods  in  Fracture  Mechanics  (Editors:  Luxmore 
and  Owen)  Proceedings  of  International  Conference  on  Numerical  Methods  in 
Fracture  Mechanics,  Swansea,  UK,  January  9-13,  1978)  pp.  52-67,  1978. 

30.  M.  Nakagaki,  W.  H.  Chen,  and  S.  N.  Atluri,  "A  Finite  Element 
Analysis  of  Stable  Crack  Growth:  Formulation  of  Criteria",  to  appear  in 
an  ASTM  STP,  Proceedings  of  ASTM  National  Symposium  on  Elastic-Plastic 
Fracture,  Atlanta,  GA,  November  1977. 

31.  S.  N.  Atluri  and  K.  Kathiresan,  "Surface  Flaws  in  Pressure  Vessels",  | 

Invited  paper.  Nuclear  Engineering  and  Design.  1978.  (to  appear).  ;■ 

32.  S.  N.  Atluri  and  K.  Kathiresan,  "Stress  Analysis  of  Typical  Flaws 
in  Aerospace  Structural  Components  Using  a Three-Dimensional  Hybrid  Element 
Procedure",  in  Proceedings  of  19th  AIAA/ASME/SAE  Structure.  Structural 
Dynamics  and  Materials  Conference.  Bethesda,  MD,  April  3-8,  1978. 
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33,  S.  N.  Atluri,  "A  Review  of  Numerical  Methods  for  Engineering 
Fracture  Analysis",  to  be  presented  and  published  at  2nd  International 
Conference  on  Fracture  Mechanics,  Bangalore,  India,  March  1979. 


5.3  PAPERS  PRESENTED  AT  PROFESSIONAL  MEETINGS 

34 . AIAA/ASME/SAE  15th  Structures,  Structural  Dynamics  and  Materials 
Conference,  Las  Vegas,  Navada,  April  1974,  paper  title  same  as  cited  in 
Ref.  1. 

35.  7th  U.S.  National  Congress  of  Applied  Mechanics,  Boulder,  Colorado, 

June  1974,  paper  title  same  as  cited  in  Ref.  2. 

36 . 1974  International  Conference  on  Finite  Element  Methods  in  Engineering, 
Cnivcrsity  of  New  South  Wales,  Sydney,  Australia,  August  29-31,  1974, 

paper  title  same  as  cited  in  Ref.  6, 

37 . NATO  Advanced  Study  Institute  on  Continuum  Mechanics  Aspects  of 
Rock  Fracture  and  Geodynamics,  Reykjavik,  Iceland,  August  1974,  presented 
an  invited  paper,  "Finite  Element  Analysis  of  Cracks  Between  Dissimilar 
Media" . 

38 . ASTM  Symposium  on  Fracture  of  High  Modulus  Fibers  and  Their 
Composites,  National  Bureau  of  Standards,  Gaithersburg,  MD,  paper  title 
same  as  cited  in  Ref.  4. 

39 . Conference  on  Fundamental  Aspects  of  Deformation  and  Fracture 

of  Composite  Materials.  Battelle  Seattle  Research  Center,  Seattle,  Washington, 
February  1975,  presented  an  invited  paper,  "Stress  Intensity  Factors  of 
Cracked  Orthotropic  Plates",  (Presented  by  A.  S.  Kobayashi) . 

40.  ASME  Joint  Western  and  National  Applied  Mechanics  Conference.  Univer- 
sity of  Hawaii,  Honolulu,  March  1975,  paper  title  same  as  cited  in  Ref.  8. 
(Presented  by  A.  S.  Kobayashi). 

41 . AICA  International  Symposium  on  Computer  Methods  for  Partial 
Differential  Equations.  Lehigh  University,  Bethlehem,  PA,  June  1975,  paper 
the  same  title  as  cited  in  Ref.  5. 

42.  AFRPL/Edwards  AFB  Contractors  Meeting,  California  Institute  of 
Technology,  Pasadena,  California,  May  1975,  invited  presentation  on  "3-D 
Cracked  Elements". 

43 . ASTM  Committee  E-24  Meeting.  9th  U.S.  National  Symposium  on 
Fracture  Mechanics.  Pittsburgh,  PA,  August  1975,  invited  presentation, 

"Analysis  of  Large-Scale  Yielding  Fracture  Problems". 

44 . Third  International  Conference  on  Structural  Mechanics  in  Reactor 
Technology,  University  of  London,  September  1975,  paper  title  same  as  cited 
in  Ref.  7. 
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45.-46.  10th  Annual  Meeting  of  Society  of  EnRineerinR  Science.  University 
of  Texas,  Austin,  Texas,  October  1975,  paper  titles  same  as  cited  In  Ref.  9 
and  10. 

47.  "Workshop  on  3-D  Fracture  Analysis",  Battelle  Columbus  Lab., 

Columbus,  Ohio,  April  1976,  invited  presentation  on  "3-D  Hybrid  Cracked 
Elements" . 

48 . 8th  South  Eastern  Conference  on  Theoretical  and  Applied  Mechanics, 

VPI  and  SU,  Blacksburg,  VA,  April  1976,  paper  title  same  as  cited  in  Ref. 

11. 

49.  17th  AIAA/ASME  SDM  Conference,  Valley  Forge,  PA,  May  1976,  paper 
title  same  as  cited  in  Ref.  12. 

50.  10th  U.S.  National  Symposium  on  Fracture,  Philadelphia,  August 
1976,  title  same  as  cited  in  Ref.  13. 

51.  13th  Annual  Meeting  of  Society  of  Engineering  Science.  NASA-Langley, 
Hampton,  VA,  November  1976.  (Ref.  16). 

52.  18th  AIAA/ASME/SAE  SDM  Conference.  San  Diego,  California,  March  1977, 
title  same  as  cited  in  Ref.  22,  (presented  by  J,  S.  Cheng). 

53.  "Crack  Elements  for  Bi-Material  Fracture,"  JANNAF  Propulsion 
Conference,  Edwards  Air  Force  Base,  CA,  March  1977. 

54.  "Inner  Surface  Flaws  in  a Thick  Walled  Pressure  Vessel,"  3rd 
International  Conference  on  Pressure  Vessel  Technology,  ASME/JSME,  Tokyo, 

Japan,  April  1977. 

55.  "Fracture  Initiation  in  Plane  Ductile  Fracture  Problems",  3rd  Inter- 
national Conference  on  Pressure  Vessel  Technology,  ASME/DSME,  Tokyo,  Japan, 
April  1977. 

56.  "Fracture  Analysis  Under  Mixed  Mode  Conditions",  2rid  ASCE  Engineering 
Mechanics  Specialty  Conference,  Raleigh,  NC,  May  1977. 

57.  "Stress  Analysis  of  Cracks  in  Elasto-Plastic  Range",  4th  International 
Conference  on  Fracture,  University  of  Waterloo,  Canada,  June  1977. 

58.  "Outer  Surface  Flaws  in  Pressure  Vessels,"  4th  International 
Conference  on  Structural  Mechanics  in  Reactor  Technology,  San  Francisco, 

CA,  August  1977. 

59.  "Analysis  of  Stable  Crack  Growth  in  Ductile  Materials",  9th  SAMPE 
National  Conference,  Atlanta,  GA  October  1977. 

60.  "Criteria  for  Stable  Crack  Growth",  ASTM  National  Symposium  on 
Elastic-Plastic  Fracture,  Atlanta,  GA  November  1977. 

61.  "Surface  Flaws  in  Plates",  14th  SES  Conference,  Lehigh  University, 
Bethlehem,  PA  November  1977. 
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62.  "Stress  Analysis  of  Typical  Flaws  in  Aerospace  Structural  Components", 
19th  AIAA/ASME,  Structures,  Structural  Dynamics  Conference,  Bethsda,  MD 

April  1978. 

63.  "Finite  Element  Analysis  of  Elastic-Plastic  3-D  Cracks",  Joint  AFOSR/ 
AFRPL  Rocket  Propulsion  Research  Meeting,  Edwards  Air  Force  Bse,  April  1978. 


5.4  SEMINARS  AT  UNIVERSITIES 

1.  "On  Elastic-Plastic  Fracture  Mechanics",  at  Department  of  Aeronautical 
Engineering,  Indian  Institute  of  Science,  Bangalore,  India,  August  21,  1974. 

2.  "Computational  Methods  in  Fracture  Mechanics",  National  Aeronautical 
Lab,  Bangalore,  India,  August  22,  1974. 

3.  "Computational  Fracture  Mechanics",  College  of  Engineering,  Boston 
University,  Boston,  Mass.,  January  1976. 

4.  "Hybrid  Finite  Elements  for  Fracture  Analysis",  four  lectures  at 
Short  Course  on  Advances  in  Finite  Element  Methods,  University  of  Tennessee 
Space  Institute,  Tullahoma,  TN  March  1978. 

5.  "3-D  Crack  Analysis  by  Finite  Elements",  University  College,  Cork, 
Ireland,  January  14,  1978. 
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FIG.  1 (Taken  From  Ref.  7) : Typical  Result  for  a center  cracked  orthotropic 
sheet;  crack  symmetric  w.r.t  axes  of  orthotropy. 
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3 (From  Ref.  6):  A Typical  convergence  study. 
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FIG.  45 : Corner  Cracks  Emanating  from 

Holes  in  Plates 


FIG,  5 (From  Refs.  2 and  21):  A Typical 


3-D  Crack  Problem. 
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Crack  Front 


No  of  Elem  = I56 


No  of  DOF* 2670 


FIG.  46:  Finite  Element  Breakdown  of 

Quarter  of  the  Corner  Cracked 
Hole  Problem 


FIG.  6 (From  Ref.  21  and  32):  Typical  Finite  Element  Mesh  for  Problem  in 


FIG.  51:  Variation  of  Stress  Intensity  Factor 
for  Corner  Cracks  Emanating  from 

Hole 

(a/c-1.5;  c /R • 1.0  ;a/t» 0.75) 

FIG.  7 (Taken  From  Ref.  32):  Typical  Result  for  a corner  crack  near  a 
Fastener  Hole. 
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SECTION  - AA 


FIG. 2.  Semi- Elliptical  Surface  Flaws  in 
Plates  in  Tension  and  Bending 

FIG.  8 (Taken  From  Ref.  28):  A Typical  3 crack  Problem  of  a surface  flawed 


I Variation  of  Stress  Intensity  Factors  for  a Semi- 

Elliptical  Surface  Crack  in  a Thin  Plate  (a/t:-8) 


FIG.  9 (Taken  from  Ref,  33':  A Typical  3-D  Results  for  surface-flawed  Tension 
Plate:  Comparison  of  Hybrid  Element  Results  with  an  alternate  Finite 
Element  Result. 


FIG  2:  Finite- Element  Model  of  3-Point  Bend  Bar 
(shown  in  insert)  and  J- Integral  Paths 

FIG.  10  (Taken  From  Ref.  IGl : A Typical  Specimen;  Analysis  of  Large-Scale 
Yielding  Fracture. 

30 


\\ 


§i 


\\\  : 

\\\  I 

'•V 

^ V 


X ^ 

Sg  < 

u.  9 

sS  t 

-J  Ul  — 


13  Si  o 

<90 

ass 

O)  It.  CD 


Lastlc  Yielding  Specimen. 


12  (Taken  from  Ref.  16):  Tjrpical  Results  for  J vs.  Gauge  Point  Displace- 
ment; Comparison  of  Various  Approximations  with  experimental  data. 


CORRESPONDING  COD  VALUES 
OBTAINED  BY  ! 

LINEAR  extrapolation 
V.'ELLS'  FORMULA 

wk-  CODA  FORMULA 

-o-  ROBINSON  ET  AL'S  FORMULA 

J - (Ty'  COD 


CRACK  TIP  OPENING 


(X  10  IN) 


FIG  14  J VS  COD  Relation  for  a 3 -Point  Bend  Bar 

FIG.  14  (Taken  From  Ref.  ISt : Correlation  of  Different  Fracture-Initiation- 
Characterizing  Parameters. 


FIG  2.  FINITE  ELEMENT  IDEALIZATION  OF 
CENTER  CRACKED  SQUARE  PLATE 
UNDER  UNIAXIAL  TENSION. 


FIG.  16  (Taken  from  Ref.  30): 
Crack  Growth. 


A Typical  Model  Problem  in  Analysis  of  Stable 
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Gfl/G 


BY  KFOURI  AND  HiLLER 


fla/ao=  .05 


FI6  3.  CRACK  SEPARAriON  ENERGT  RELEASE  RATES  : 
CONSTANT  aa  AT  VARIOUS  LOAD  LEVELS. 


IG,  17  (Taken  From  Ref.  29'':  Study  of  crack  Separation  Energy  Release  Rates 
'n  Ductile  Materials:  The  parameter  S is  the  ratio  of  Growth  Stop 
to  the  Plastic  Yield  zone  size. 


CONSTANT  40  AT  VARIOUS  LOAD  LEVELS. 


FIG.  18  (Taken  From  Ref,  30"':  Further  Elaboration  of  Results  in  Fig.  17. 


CRACK  GROWTH  CURVE. 


P/2B  (MN/m)  (KPom)  J (KPam) 


FIG  7.  VARIATION  OF  J,  G*^AND  P/2B 
DURING  CRACK  GROWTH. 


FIG.  20  (Taken  From  Ref.  30):  Results  for  Various  Physical  Parameters  that 
may  govern  stability  of  Crack  Growth,  in  the  Simulation  of  Data  in 
Fig.  19. 
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ENERGY  RATE  (KPom) 


.2  .4  .6  .8  10  1.2 


Q- Qo  (mm ) 

FIG  9.  VARIATION  OF  G)^,6'^.AND  ENERGY 
DISSIPATION  IN  PROCES  ZONE  DURING 
CRACK  GROWTH. 

FIG.  21  (Taken  Fig.  Ref,  30);  Continuation  of  Results  Similar  to  those  in 
Fig.  19. 


41 


CTOA  (Rad) 


I 


.2  .4  .6  .8  1.0  12 

a-  do  ( mm  ) 

FIG  10.  VARIATION  OF  CRACK -TIP  OPENING 
ANGLE  DURING  CRACK  GROWTH. 

FIG.  22  (Taken  From  Ref.  30):  Continuation  of  Results  Similar  to  those  in 
Fig.  19. 
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